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Abstract 

We discuss predictions for the interaction energy of the fundamental monopoles 
in gluodynamics introduced via the 't Hooft loop. At short distances, the heavy 
monopole potential is calculable from first principles. At larger distances, we ap- 
ply the Abelian dominance models. We discuss the measurements which would 
be crucial to distinguish between various models. Non-zero temperatures are also 
Q-i! considered. Our predictions are in qualitative agreement with the existing lattice 

O ' data. 



1 Definition of the heavy monopole potential 



In this note we discuss the interaction of the fundamental monopoles in gluodynamics. 
The fundamental monopoles can be introduced via the 't Hooft loop |TJ , and are best un- 
derstood on the lattice. The monopoles correspond to point-like objects in the continuum 
limit and can be visualized as end-points of the Dirac strings which in turn are denned 
as piercing negative plaquettes. In more detail, consider the standard Wilson action of 
SU(2) lattice gauge theory (LGT): 

Si at (u) = -pYl h±u p , (1) 

where f3 = 4/g 2 , g is the bare coupling, the sum is taken over all elementary plaquettes 
p and Up is the ordered product of link matrices in the fundamental representation along 
the boundary of p. Then the 't Hooft loop is formulated (see, e.g. 0, U and references 
therein) in terms of a modified action S(f3, —/3): 

S(p,-0) = -P Y, ! Ti[/ p + ^ l^U pj (2) 
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where *E 3 - is a manifold which is dual to a surface spanned on the monopole world-line j. 
Introducing the corresponding partition function, Z(/3, —(3) and considering a time-like 
planar rectangular T x R, T ^> R contour j one can define 

V mrh (R) = -f^^jy- (3) 

Since the external monopoles become point-like particles in the continuum limit the 
potential V mm (R) is the same fundamental quantity as, say, the heavy-quark potential Vqq 
related to the expectation value of the Wilson loop. By analogy, we will call the quantity 
V mifl {R) the heavy monopole potential. Because of the asymptotic freedom, one would 
expect that at short distances the potential V m m(R) is predictable from first principles. 
Some preliminary work is needed, however, since originally V m m(R) is formulated in the 
lattice terms. The continuum version of the 't Hooft loop was worked out in [Q . Moreover, 
it turns possible || to reformulate the problem in the Lagrangian approach a la Zwanziger 
U . The crucial point is that the role of the Dirac string, explicit in the lattice construction 
(0), reduces to particular rules of the ultraviolet regularization. The papers @, |5| set up 
a theoretical framework which we will exploit here. 

Recently, first direct measurements of V mm ,(R) on the lattice were reported [fj, |3j. The 
main conclusion is that at large distances the potential V m m(R) is of the Yukawa type. 
Motivated by these measurements we explore in this note how far one can reach with 
predicting the heavy monopole potential theoretically. As is already mentioned above, we 
are not trying to develop new models but rather to elucidate the relation between various 
theoretical hypotheses and measurable properties of the potential V m m{R)- 



2 Potential at zero temperature 
2.1 Coulomb-like potential 

Consider first the potential V m m{R) at short distances. Then the interaction between the 
monopoles is Coulomb-like (see 0, || and references therein): 

(27T \ 2 1 7T 1 

J(R)) A^R = ~J{K)R' (4) 

where g(R) is the running coupling of the gluodynamics. Eq. (f|) makes manifest that 
monopoles in gluodynamics unify Abelian and non-Abelian features. Namely, the overall 
coefficient, ir/g 2 is the same as in the Abelian theory with fundamental electric charge 
e = g while the running of the coupling g 2 is governed by the non-Abelian interactions. 

Prediction (|J) is rooted in the very foundations of the theory. Namely, the U(l) 
overall normalization is based on the classification of the monopoles in non-Abelian gauge 
theories (see, e.g., [0] and references therein). The running of the coupling reflects the 
general rule that the effect of the fluctuations at short distances can be absorbed into the 
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renormalization of the coupling. What might be not so evident, is which distances can be 
considered small in practice. We shall come to discuss this point later. 

To derive ([|) explicitly one establishes first the classical equations of motion corre- 
sponding to the 't Hooft loop: 

d 2 A^ = -— ^(n^), (5) 

where is the world sheet of the Dirac string whose end points represent the monopoles, 
star symbol means the duality operation and n a is a unit vector in the color spaceQ whose 
direction is a matter of gauge fixing. In particular, there exist gauges in which n a is 
constant. 

Eq. d5|) demonstrates that only the boundary of the Dirac string determines the in- 
teraction. Therefore, one can hope to introduce Lagrangian which would reproduce the 
same results in terms of particles alone, without strings. A well known example of such a 
construction is the Zwanziger Lagrangian || which allows to describe photons interacting 
with both magnetic and electric charges. In case of gluodynamics the Zwanziger-type 
Lagrangian was found in j5|: 

L dnal (A a ,B) = - A {F* v f + ±[m-(dAB-z*n a F a )} 2 + ZjjB, (6) 

where j is the magnetic current, is the non-Abelian field strength tensor and m M is 
an arbitrary constant vector, m 2 = 1. Note that the vector m enters also the original 
Zwanziger Lagrangian in the QED case. The vector field B can be called dual gluon, 
although because the A a and B fields are mixed the number of the physical degrees 
of freedom corresponding to (^) is the same as in the gluodynamics itself. One can 
demonstrate that the Lagrangian @ reproduces the heavy monopole potential @). The 
Dirac string, however, re-emerges in a way, through rules of the ultraviolet regularization 
which uniquely define the theory (^) on the quantum level. 

A salient feature of the Eq. (Q) is that the dual gluon B is a U{1) particle despite of 
the fact that the "ordinary" gluon are in an adjoint representation of the non-Abelian 
group. There is no violation, however, of the color symmetry through the mixing of the 
B and A a fields since the choice of n a is a matter of gauge fixing, see for details j|, [5] . In 
the simplest form, the argument is that an averaging over all possible embeddings of the 
£7(1) into the non-Abelian group is understood. Thus, by checking (|^) one would verify 
that the dual gluon in gluodynamics is indeed a U{1) gauge boson. 

The message so far is that at one-loop level the interaction of the heavy monopoles 
is understood theoretically no worse than the interaction of the heavy quarks. There is 
a word of caution, however, concerning higher loops. The point is that the monopole 
field is proportional to g~ x . As a result, any number of interaction of virtual gluons with 
the external monopole field is equally important and the evaluation of the corresponding 
determinant looks a very difficult problem, with no analytical solution at all ||. The 

1 For simplicity we consider only the SU (2) gauge group. 
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leading log approximation is an exception to the rule f| and is given by the simple graph 
with insertion of two external fields. Thus, it may take long before the next term in the 
expansion in g 2 (R) is found. 



2.2 Power corrections 

As far as the physics of short distances is concerned, the next step is to consider power 
corrections to (^). Theoretically, the prediction is that there is no linear in R term at 
small R: 

71 1 

V mm (R) « -g2(fy R + a o A QCD + a 2 A 3 QCD R 2 + ... (7) 

where ao,2 are constants sensitive to the physics in the infrared. To the contrary, the 
coefficients in front of the odd powers of R are sensitive to the short distances. The same 
is true in case of the heavy quark potential and we refer the reader to || for the reasoning 
and further references. The absence of the linear term follows then from the absence of a 
dimension d = 2 quantity in the gluodynamics. 

Note that in case of the heavy quark potential this logic can in fact be challenged 
[ [Toll . Indeed, the power corrections are obviously sensitive to the vacuum structure. The 
picture that we have in mind is that monopoles condense in the vacuum (see below). 
Then the external quarks are connected by a Dirac string and since the string is infinitely 



thin, the physics in the ultraviolet can change lOfl . The case that we are considering now, 
that is, magnetically charged external probes embedded into the vacuum with monopoles 
condensed, is much simpler. Indeed, in the dual language we do not need then Dirac 
strings at all. Thus, there is no mechanism to change (Jfty. Things would change if we 
included light quarks. Then the quark condensate could invalidate (|7]) since the quarks are 
carrying color. However, the dynamical quarks would also make visible the Dirac string 
connecting the external monopoles, and the whole construction would radically change. 
We would not consider that case in detail. 

To summarize, the absence of the linear in R term at short distances follows also 
from the first principles and eventually is related to the fact that the monopoles are not 
confined. 



2.3 Yukawa potential 

It is a common point that at large distances V mm ,(R) is of the Yukawa type, see, e.g., 
[0,0,1,0,0: 

V m m(R) = — C ■ — — . (8) 

a 

In fact, it is a manifestation of the general statement that either electric or magnetic 
charges are confined JOJ. At a closer look, however, there is a variety of model depen- 
dent predictions for the parameters C and [l. Thus, experimental determination of these 
parameters could distinguish between various models and we will discuss the models one 
by one: 
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(i) Continuing with our discussion of the fundamental gluodynamics, Eq. (0) tells us 
that the Yukawa potential (§) cannot be true at all the distances. Indeed, absence of the 
linear term is inconsistent with (H). In other words, there is no reason to believe that the 
Yukawa potential matches smoothly the Coulomb-like potential and, generally speaking, 
C/ir/g 2 . Thus, there are no model independent predictions either for ii or C. 

(ii) Historically, the first prediction for fi was obtained in ||11|| . Namely, it was shown 



that to any order in the strong coupling expansion the mass fx coincides with the mass of 
0++ glueball: 

/i = m G . (9) 

As for the constant C, there is again no reason for it to be the same for the Yukawa and 
Coulomb like potentials so that, generally speaking, C ^ ^ j g 1 . 

(Hi) The Yukawa potential fl8|) at all the distances arises naturally within effective field 
theories with monopole condensation^. To describe the condensation in the field theoret- 
ical language, one introduces a new (effective) field <fi which is interacting minimally with 
the dual gluon (for review and references see, e.g., || 14|). Consider first the Lagrangian 
proposed in 0. 

Leff = L dual (A a , B) + S mggs (B , <f) m ) . (10) 

where Lduai is given by Eq. (§) and -^Higgs is the standard Higgs part of the Abelian Higgs 
model action. The vacuum expectation value of the Higgs, or monopole field is, of course, 
of order Aqcd- Within this model, 

7T 

H = m v , C = — , (11) 

g 2 

where my is the mass of the vector field B acquired through the Higgs mechanism. An 
attractive feature of the model ( |T0"D is that it can incorporate the Casimir scaling ||. 



It is worth emphasizing that the predictions (|TT| ) differ substantially from the cases 
(i), (ii) discussed above. First, my may not coincide with any glueball mass. This is not 
in contradiction with the general principles since there is no spectral representation for 
the correlator of two magnetic currents j (which are sources of the field B) [|J. Indeed, 
although the currents do not carry any color index explicitly they are not color singlets 
either since their color orientation is a matter of a gauge fixing. Moreover the monopole 
trajectory j is to be understood as the boundary of the Dirac-string world sheet. The 
Dirac string is infinitely heavy in the continuum limit so that all the intermediate states 
are in fact infinitely heavy. 



As for the magnitude of my, the estimates usually give my ~ 1 GeV, see, e.g., [ [0| . 
More generally, if it turns out that fi is indeed smaller than the lowest glueball mass, it 
would be a serious argument in favor of the Abelian dominance models. If the value of 
fx does not distinguish between the models, then it would be crucial to check that the 
potential is vector like, as is predicted by the model (|10|). 

2 Note that the monopoles which condense are of course not the fundamental monopoles which are 
introduced via the 't Hooft loop as external probes. Instead, the monopoles "living" in the QCD vacuum 
have a double magnetic charge. 
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At very short distances the potential @ should yield to the potential (0) obtained 
within the fundamental gluodynamics. However, there exist various pieces of evidence 
that the effective theory ( fLOf) is valid down to such small scales that in practice it covers 
all the distances available for the lattice measurements nowadays [|K], [L6|. Clearly, a 
careful study of the V mm (R) could be crucial to confirm or reject these speculations. 

(iv) It is worth emphasizing that the value C = n/g 2 is specific for the model (|T0|) 
and is not true in a generic Abelian dominance model. To use an analogy, consider the 
world with spontaneously broken chiral symmetry of strong interactions and with weakly 
interacting vector bosons added. Then the model fllPl) would be analogous to assuming 
that the massless pions are coupled to the W bosons but do not interact directly with 
the nucleons. Thus, we could have added to ( |10|) interaction of the Higgs bosons with the 
fundamental monopoles. In other words, since the monopoles are condensed the monopole 
charge of particles is not well defined and the constant C is not constrained, generally 
speaking. 



3 Non-vanishing temperatures 
3.1 Dimensional reduction approach 

The screening mechanism at high temperatures is the Debye screening. As is noted above 
the classical limit of the state created by 't Hooft loop is an Abelian monopole pair. Thus 
it is natural to evaluate the Debye mass within the Abelian dominance models. Note that 
even if one assumes the Abelian dominance, it is an open question which U(l) is to be 
used to classify the monopoles. At short distances, there is no such problem since the 
coupling g 2 is small and the only subtle point is the averaging over all the embeddings 
of the U(l) into 577(2). At large distances when g 2 ~ 1, different definitions of the U(l) 
result in different physics. We assume that the maximal Abelian gauge is used. 

Our estimates of the mass p{T) M at high T are based on the observation that 
the Abelian model which corresponds to the high temperature gluodynamics is the 3D 
compact U(l) theory. Therefore, at high temperatures the screening mass p coincides 
with the corresponding Debye mass: 

/j, 2 = m 2 D = 16tt 2 4, (12) 

e 3 

where p is the density of monopoles and e% is the corresponding three-dimensional coupling 
constant. To estimate the temperature dependence of m D we use the numerical results of 
Ref. [17|, where the density of Abelian monopoles was obtained^: 



p = 2~ 7 (1 ±0.02)e°. (13) 



Note that the original result of Ref. ||17[| refers to the lattice monopole density pi at — 0.50(1) 0q, 



where (3q is a three dimensional coupling constant which is expressed in terms of the 3D electric charge 



e3 and lattice spacing a as 0q = 4/(ae§). The physical density p of monopoles is given by p = pi a t. a 
which can easily be transformed into Eq. (|l3|). 
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At high temperatures we can use the dimensional reduction formalism and express 3D 
coupling constant in terms of 4D Yang-Mills coupling g. At the tree level one has 



el(T)=g 2 (A,T)T, (14) 
where g(A,T) is the running coupling calculated at the scale T, 

S " (A ' T) = IsWx) + 4iW 21 °s(x)] ■ (15) 

and A is a dimensional constant which can be determined from lattice simulations. 

At present, the lattice measurements of the A parameter are ambiguous and depend 
on the quantity which is used to determine it. We use two "extreme" values of A. Namely 
in Ref. |0| the lattice data for the gluon propagator have been used to determine the so- 
called "magnetic mass" in high temperature SU(2) gluodynamics. These measurements 
indirectly provide, 

A = 0.262(18) T c , (16) 
where T c is the temperature of the deconfmement phase transition, T c pa 0.69-y/cx (see, 



e.g., fll9|). In Ref. P0| , on the other hand, a the spatial string tension has been calculated 



and a smaller value of A was found: 

A = 0.076(13) T c . (17) 

Another uncertainty in our prediction is that we have used the dimensional reduction 
which is supposed to work well only at asymptotically high temperatures. In practice, 
however, the temperatures used in lattice measurements are only few times larger T c . 

3.2 Monopole potential in Abelian projection 

So far we assumed that the heavy monopole potential is measured via the 't Hooft loop 
in full, non-Abelian lattice simulations. However, within the Abelian dominance model 
all the interactions are described by QED-like interactions of the magnetic and electric 
charges. Moreover, the magnetic currents in the vacuum can be measured directly. Then 
the potential energy of external sources can, in principle, be evaluated using the ensembles 
of the magnetic currents. In case of external color charges, or Wilson loop such an 
approach is well known. Also, the intermonopole potential has been studied on the lattice 
both at zero |yj and non-zero |12| temperatures. In this section we address this issue 
anew and emphasize that in fact calculation of the heavy monopole potential involves 
extra parameters. However, the situation is simplified greatly in the limits of very low 
and very high temperatures. 

The Abelian monopole action in the Maximal Abelian projection depends on many 
parameters. The action contains the term 5 ,< - c ' ou ^ = const. A~ 1 *j), responsible 
for the Coulomb exchange between monopoles as well as various n— point, n > 2, local 



monopole interaction terms, Ref. |24||. Here A 1 is the lattice inverse Laplacian, *j is 



7 



the monopole current on the dual lattice. Throughout this Section we are using the 
differential form notations on the lattice p5[ 



Let us consider the following simple monopole action with 2-point (self) interaction 
as an example: 

S ^™)[*j] = 47r 2 /?(*j, A~ u j) +47r 2 7|rj|| 2 , (18) 

This action corresponds to the London limit of the (dual) lattice Abelian Higgs model |26f 
in which the role of the Higgs field is played by the monopole field. Another representation 
of model fll8"D is the compact gauge field representation p4j : 

7T 

Z = J V6e- s(comp, ( d9 \ 
e- s(mP> W = exp{-/3||d^ + 27™|| 2 - 7 (d0 + 2™, A(d0 + 2™))} , (19) 

where 9 is the compact gauge field and n is the integer valued auxiliary plaquette field. 
The intermonopole potential can be studied then with the help of the Abelian 't Hooft 
loop, 

Hf j q = exp{S {comp) (d9) - S icomp \d9 + 2nq *Sj)} , (20) 

where the surface £j ends on the trajectory of an external monopole with charge q. 

Similarly to eq. (|3|) the quantum average of the 't Hooft operator gives the static 
monopole potential: 



mon 



yab ( R ) = -i ln§— , Z™ on = E e- 5(m ° n) ^* J ], (21) 

mm, q \ ) rj-i nmon ' 9 i. < \ ' 



T Zl 



8*j=0 



where the summation is taken over all closed monopole trajectories. The 2-point inter- 
action term provides the local interaction between external J and dynamical j currents, 
Sint(*j, *J) = 8g7r 2 7(*j, * J). In our case the external current consists of two disconnected 
pieces, *J = *J m + *J* m , separated by the distance R. 

The term Si nt (j, J) is non-zero iff the dynamical monopole current touches the ex- 
ternal monopole current J, see Figure |](a). Thus, the 2-point interaction terms affect 
the intermonopole potential provided the dynamical monopoles overlap locally with the 
external sources. 

Let us discuss the role of the local n-point interaction terms at the finite temperature. 
At sufficiently small temperatures the system is in the confinement phase in which the 
entropy of the monopole currents prevails over their energy. Thus there exist arbitrarily 
long dynamical monopole currents j connecting the external (anti-) monopole trajectories 
J m (J m ), Figure 0(a). In other words, the monopole currents in the confinement phase of 



gluodynamics are percolating |27[]: the probability for two distinct points to be connected 



by a monopole trajectory does not depend on the distances between them at sufficiently 
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large separations. Since the dynamical currents are not getting suppressed as their length 
is increased the distance dependent part of the inter-monopole potential is not affected 
by the n-point interaction terms. 




Figure 1: The static external currents J m and J m and the dynamical monopole currents 
j at confinement, (a), and the deconfinement, (b), phases. 



As temperature increases the system goes into the deconfinement phase. Now the 
n-point terms become important since the monopole currents acquire non-zero free en- 
ergy per unit monopole trajectory length. However at sufficiently high temperatures the 
dynamical monopole trajectories become static (i.e., wrapped on the compactified direc- 
tion with a minimal length), see Fig. [l|(b), and n-point terms do not affect the distance 
dependence of the correlators for well separated external monopoles. 

To summarize we expect that the magnetic screening mass obtained using only the 
long range part of the action should be correct at asymptotically small and large temper- 
atures while at the intermediate temperatures the result is modified by the local current 
interactions. 

3.3 Comparison with the data 

Data on V mm (R) are only beginning to accumulate and we summarize briefly the com- 
parison of the lattice measurements with predictions above: 

Short distances. The Coulomb-like potential @ is confirmed in the numerical simu- 
lations in the classical approximation |J. There is no running of g 2 on this level of 
course. As for the full quantum simulations the normalization (|]) of the potential at short 
distances is confirmed within a factor of about 2, Ref. [f2"5|]. 
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As for the power corrections, all the data so far |2|, ||| are fitted smoothly with a 
Yukawa potential @. Thus, there is no evidence whatsoever that the linear term at short 
distances is to be subtracted. However, no explicit bound on the linear term at short 
distances has been obtained either. 

Screening mass at T = 0. Existing data @, |J seem to be not accurate enough to 
distinguish between the models (|9p, ([XT]) . Moreover, no checks of the vector-exchange 
have been made. 

Temperature dependence of the screening mass. On Figure we have summarized 
the existing data on temperature dependence of the screening mass together with our 
predictions, Eqs. (0-jl5)- The direct measurement of the screening mass [§ is shown 
by the diamonds and squares, respectively. The character of the predicted temperature 
dependence is reproduced by the data. Although numerically the screening mass seems to 
be systematically higher than is predicted, the prediction based on the Abelian dominance 
hypothesis is rather close to the numerical results within the errors. 




Figure 2: Temperature dependence of the screening mass in the heavy monopole potential 
(H). Shaded region denotes the theoretical prediction, Eqs. (|ll|), (|12| - [T5| ). Symbols denote: 
x - the dual gluon mass my « lGeV as reported in Ref. [Tj|; * - the mass of the lightest 
++ glueball |1|; y - magnetic screening mass measured in Maximal Abelian gauge, 



Ref. [12]; O and □ - the screening masses measured at various temperatures in Ref. 



and Ref. respectively. 
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The triangles on Figure 1 denote the results of Ref. [|12|] obtained in the Maximal 
Abelian projection using the method described in the previous Section. At small tem- 
peratures these results are quantitatively in agreement with Refs. [[| |3| while at higher 
temperatures the screening mass obtained in Ref. |12| falls essentially lower than is indi- 
cated by the measurements in the full gluodynamics @, §] . It is worth emphasizing that 
the calculations [0 took into account only the Coulomb-like interaction between the 
monopoles. As is explained in the preceding section, this is justified at very low and very 
high temperatures, but not at the intermediate temperatures. This observation allows 
to understand why the results of 0, |3| and |L2| agree at T = and tend to disagree at 
intermediate temperatures. 



Conclusion 

To summarize, the heavy monopole potential at small distances is of the Coulomb type 
with a known overall normalization. The large-distance potential is of the Yukawa type. 
The descent of the interaction at large distances and at non-zero temperatures is quali- 
tatively consistent with the monopole dominance model. Thus the monopole dominance 
model of SU(2) gluodynamics based on the dual superconductivity predicts the inter 
monopole potential which is consistent with the numerical data both at small and large 
distances at high as well as at low temperatures. Further data are desired to allow for a 
detailed quantitative comparison. 
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